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Abstract 



. The paper deals with a spatially homogeneous and totally anisotropic Bianchi II cosmological models 

representing massive strings in normal gauge for Lyra's manifold. The modified Einstein's field equations 
have been solved by applying variation law for Hubble's parameter. This law generates two type of solutions 
for average scale factor, one is of power law type and other is of exponential law type. The power law 
describes the dynamics of Universe from big bang to present epoch while exponential law seems reasonable 
to project dynamics of future Universe. It has been found that the displacement actor (/?) is a decreasing 
function of time and it approaches to small positive value at late time, which is collaborated with Halford 
(1970) as well as recent observations of SN la. The study reveals that massive strings dominate in early 
Universe and eventually disappear from Universe for sufficiently large time, which is in agreement with the 
current astronomical observations. 
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1 Introduction 

>' 

In recent years, there has been considerable interest in string cosmology because cosmic strings play an important 
role in the study of early Universe. It is generally assumed that after the big bang, the Universe may have 
undergone a series of phase transitions as its temperature lowered down below some critical temperature as 
predicted by grand unified theories [_] — [_!]. At the very early stages of evolution of the Universe, it is believed that 
during phase transition the symmetry of the universe is broken spontaneously. It can give rise to topologically 
stable defects such as domain walls, strings and monopoles. In particular, cosmic strings are produced in the 
breaking of U(l) symmetry, are good candidates to seed the formation of galaxies. These cosmic strings have 
stress-energy, and couple to the gravitational field. Therefore, it is interesting to study the gravitational effects 
. ■ that arise from strings. The pioneering work in the formulation of the energy-momentum tensor for classical 
massive strings was done by Letelier [7] who considered the massive strings to be formed by geometric strings 
with particle attached along its extension. Letelier [8] first used this idea in obtaining cosmological solutions in 
Bianchi-I and Kantowski-Sachs space-times. Stachel [9] has studied massive string. 

The observed Universe is satisfactorily described by homogeneous and isotropic models given by the FRW 
space-time. But at smaller scales, the Universe is neither homogeneous and isotropic nor do we expect the 
Universe in its early stages to have these properties. The anomalies found in the cosmic microwave background 
(CMB) and large scale structure observations stimulated a growing interest in anisotropic cosmological model of 
Universe. Here we confine ourselves to models of Bianchi-type II. Bianchi type-II space-time has a fundamental 
role in constructing cosmological models suitable for describing the early stages of evolution of Universe. Asseo 
and Sol [10] emphasized the importance of Bianchi type-II Universe. 

Roy and Banerjee [11] have dealt with locally rotationally symmetric (LRS) cosmological models of Bianchi 
type-II representing clouds of geometrical as well as massive strings. Wang [12] studied the Letelier model in the 
context of LRS Bianchi type-II space-time. Recently, Pradhan et al. p__| __4] and Amirhashchi and Zainuddin 
[15] obtained LRS Bianchi type II cosmological models with perfect fluid distribution of matter and string dust 
respectively. Belinchon [161 117] studied Bianchi type-II space-time in connection with massive cosmic string 
and perfect fluid models with time varying constants under the self-similarity approach respectively. Recently, 
Kumar [IB] has investigated string cosmological models in Bianchi type-II space-time. 
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In last few decades there has been considerable interest in alternative theories of gravitation. The most 
important among them being scalar-tensor theories proposed by Lyra [TjJ . Lyra suggested a modification of 
Riemannian geometry which may also be considered as a modification of Wey's geometry. In Lyra's geometry, 
Weyl's concept of gauge, which is essentially a metrical concept, is modified by the introduction of a gauge func- 
tion into the structure less manifold. This alternating theory is of interest because it produces effects similar to 
those produced in Einstein's theory. Also vector field in this theory plays similar role to cosmological constant in 
general relativity [5UJ HI] ■ Several authors Sen and Vanstone [55] , Bhamra [53] , Singh and Singh [53] , Rahaman 
et al [25l [26] , Pradhan ct al [27] — [31] , Yadav [32] and Yadav et al [33] [34] have studied cosmological models 
based on Lyra's manifold in various physical contexts. 

Recently, Agwarwal, Pandey and Pradhan |35j have studied Bianchi II string cosmological model in Lyra's 
manifold. They have investigated time varying vector field by assuming direction of string along x-axis. In this 
paper, we have assumed the direction of string along z-axis, collaborated with Saha [36) . The paper is organized 
as follows. The metric and the field equations are presented in Section 2. Section 3 deals with exact solutions of 
the field equations with strings fluid. Physical behavior of the derived models are elaborated in detail. Finally, 
in Section 4, concluding remarks are given. 



2 The metric and field equations 

We consider totally anisotropic Bianchi type-II line element, given by 

ds 2 = -dt 2 + A 2 (dx - zdy) 2 + B 2 dy 2 + C 2 dz 2 , (1) 

where the metric potentials A, B and C are functions of t alone. This ensures that the model is spatially 
homogeneous. 

The field equations ( in gravitational units c = l,8nG = 1), in normal gauge for Lyra's manifold, obtained 
by Sen [37J as 

-glR+-^ --« 

where is the is the displacement field vector defined as 



4 - -zdR + - -jM k = -t;, (2) 



& = (0,0,0, /?(*)) (3) 
and the other symbols have their usual meaning as in Riemannian geometry. 

Einstein tensor. The energy-momentum tensor for a cloud of massive strings and perfect fluid distribution is 
taken as 

T] = + pXvj + pg) - Xx l x 3 , (4) 

where p is the isotropic pressure; p is the proper energy density for a cloud strings with particles attached to 
them; A is the string tension density; v % = (0, 0, 0, 1) is the four-velocity of the particles, and x 1 is a unit space-like 
vector representing the direction of string. The vectors v l and x l satisfy the conditions 

ViV 1 = -XiX 1 = -1, v l Xi = 0. (5) 

Choosing x % parallel to d/dz, we have 

x i = (o,o,cr\o). (6) 

Here the cosmic string has been directed along z-direction in order to satisfy the condition — T 2 2 . As a 
result, the off-diagonal component of Einstein tensor, viz., G\ = z(T 2 — T^) vanishes. A detailed analysis about 
the choice of energy- momentum tensor for Bianchi type-II models is given by Saha [36] . 

If the particle density of the configuration is denoted by p p , then 

p = p p + A. (7) 

The modified Einstein's field equations ([2]) for line element ([T]) with energy- momentum tensor ([4]), lead to 
the following set of independent differential equations: 

B C BC 3 A 2 3 2 _ 

B + C + BC 4B 2 C 2 4^ ~ Pl [ ' 
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C A CA 1 A 2 3 a _ 
C + I + CA + _ 4^ 



A B AB 1 A 2 



A B AB AB 2 C 2 4 

CM 1 A 2 3 
AB + ~BC + OA ~ 1B 2 C 2 + 4 



/3 2 



(9) 
(10) 

(11) 



Here, and in what follows, an over dot indicates ordinary differentiation with respect to t. The energy 
conservation equation T lj '- = 0, leads to the following expression: 



. , , A B C\ ,C n 



(12) 



and conservation of R. H. S. of cq. © leads to 

1 



-MR 



Equation P^|) reduces to 

6'r?' 



d<j>> 



Equation (fT5]) is identically satisfied for i = 1,2,3. for 



g(.9 44 ^4) 
<9x 4 



44 



d(t>4 
dt 



d(j) k 
dx 3 

= 4, eq. (TT5[) reduces to 
3 



^4r 44 

3 



fe 1 



6 4 r 4 



= 



4 44 i4 



at 



6 4 r 4 



44 



= 



which leads to 



3 • 3 , / A B 
2^ + 2 P [a + B 



= 



(13) 



(14) 



(15) 



(16) 



Thus eq. (TT21 combined with eq. (TTB1 is the resulting equation when energy conservation equation is satisfied 
in the given system. 



3 Solutions of the Field Equations 

Equations (l5])- ([TTj) and (|16l) are five equations in seven unknown parameters A, B, C, p, p, A and f3. Two 
additional constraints relating these parameters are required to obtain explicit solutions of the system. 

First, we utilize the special law of variation for the Hubble's parameter given by Berman |38j . which yields 
a constant value of deceleration parameter. Here, the law reads as 

H =l(ABC)~%, (17) 

where t > and n > are constants. Such type of relations have already been considered by Berman and Gomide 
[55] for solving FRW models. Later on, many authors (see Kumar and Singh [50], Pradhan et al [51], Yadav [52] 
and references therein) have studied Bianchi type models by using the special law for Hubble's parameter that 
yields constant value of deceleration parameter. 

Considering (ABC)s as the average scale factor of the anisotropic Bianchi-II space-time, the average Hubble's 
parameter may be written as 

»=\(i44\ 

Equating the right hand sides of (fl~3)) and (|14l) . and integrating, we obtain 

ABC = (ntt + a)", (n^O) (19) 
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ABC = c\e u \ (n = 0) (20) 

where c\ and C2 are constants of integration. Thus, the law (|17p provides power-law (IT91 and exponential-law 
([20]) of expansion of the Universe. 

Following Pradhan and Chouhan [?T], we assume that the component a\ of the shear tensor a° i is propor- 
tional to the expansion scalar {&). This condition leads to the following relation between the metric potentials: 



A=(BC) m ; 

where m is a positive constant. 

Now, subtracting (J9j) from (|8]), and taking integral of the resulting equation two times, we get 



B 

— = C4 exp 



1 



{AcS liC dt \ dt + C \l ABC 



dt 



(21) 



(22) 



where c 3 and C4 are constants of integration. In the following subsections, we discuss the string cosmology 
using the power-law (fT9|) and exponential-law (|20|) of expansion of the Universe. 



3.1 Lyra's Cosmology of String fluid with Power-law 

Solving the equations (|TB| . (TT2"j) and p^|) . we obtain the metric functions as 

A(t) = (nit + a)^^ , 



B(t) = c 4 (n£t + ci) exp 



C{t) = c^(n£t + ci)^ttt exp 



(m+ l) 2 
2PM 



(nit + ci)^T +2 + °^ 3) + ci) 2 



(^ ( ,, <t + cl) l««._iL_ ( „ H + cl) - 



(23) 
(24) 

(25) 



where Af = (9m — 3 + mn + n)(3m — 3 + mn + n) and n ^ 3. Thus, the metric ([l} is completely determined 
Equation (fl6|) leads to 

/3 I A B C 



(26) 



which reduces to 



Integrating equation (|2"T)) . we obtain 



where fci is the constant of integration. 



3/ 



P = 

/3 (nft + ci) 



/3 = fci (n# + ci)~ 



(27) 
(28) 



The expressions for the isotropic pressure (p), the proper energy density (p), the string tension (A) and the 
particle density (p p ) for the above model are obtained as 



3£ 2 [n(m + 1) - 3(m 2 - m + 1)] 



3(ro + l)(n+l) 



6(1-'") 



+ 



(m + l) 2 
3c 3 £(l ~ 2m 



{nit + Cl )- 2 + ' 1 ; y (nft + ct) ^nt 

4(9m — 3 + mn + n) 



m + 1 
2c 3 (m+ 1) 



(m+ l) 2 



£ 2 (9m — 3 + mn + n) 2 



.-(n^ + ci) x l -c^(n^ + ci)-« - 7fc 2 (n^ + ci)~", 

£(9m - 3 + mn + n) x A 4 v ; 



(29) 



9£ 2 m(2 — m) . „ . 9 15 — 33m — ran — n , „ . e(i-m) 
-(n^ + ci) _2 + — -(n^t + ci) 



+ 



(m + l) 2 
3c 3 £(l - 2m) 



m + 1 
2c 3 (m + 1) 



{nit + ciY 



4(9m — 3 + mn + n) 
_i (m+1) 2 



£ 2 (9m — 3 + mn + n) 2 



(n& + Ci)~ 



3 + mn + n) 



[nit + a) 



cl(n£t + ci) ~ + -k 2 (n£t + ciY 



(30) 
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A = 3^(2m-l)(3-n) + a + ^ + ^ 
m + 1 

3l 2 [3m(2-m)) + (3-n)(l-2m)(m+l)] ^ „,_ 2 3(13 - 35m - 3mn - 3n) „ ^Sfegft 

Pp - (m + 1) 2 [rut + d) + 4(9m _ 3 + mn + n) inn + dj 

3cs^(l — 2m) , s s , (m + 1) 2 , s r igg-m) „i 

+ ^ : ^ nit + ci " -- 1 - \ ^ rcft + ci -^^™- 2 J 

m+1 r(9m-3 + mn + n) /l 

2C3(m + 1} -(ntt + c^-tCT" 1 ] - c 2 (n£t + d)"* + ?fc 2 (nft + d)"*. (32) 



£(9m — 3 + mra + n) 4 
The critical energy density (p c ) and density parameter (ft) are given by 

Pc = U 2 (n£t + ci)" 2 (33) 

_ 3m(2 — m) 15 — 33m — mn — n . „ . gq-m) . 2 

° = 7^TlF + 12(9m-3 + mn + nK 2 ( ^ + Cl) ^ 



c 3 £(l-2m), . , 3,, (m + 1) 2 , „ , ma-m) 

+ ^ + IT (nft + ~ 3^(9m - 3 + 1 + n) 2 ^ + ^ ^ 



-4 



2C3(m + 1} -(nft + Cl )-t™- 3 ] - i-JnU + Cl )- +2 + £(»ft + c^ 2 , (34) 



3l 3 (9m-3 + mn + n) v 3£ 2 y ' 4£ 2 

The average Hubble's parameter, expansion scalar, shear scalar spatial volume (V) and anisotropy parameter 
(A)spatial volume (V) and anisotropy parameter (A) of the model are, respectively given by 

H=^(H x + H v + H z )=£(n£t + c 1 )- 1 , (35) 

6 = 3H = 3£(n£t + ci)" 1 , (36) 

o 2 = ±[(H x -H y ) 2 + (H y -H z ) 2 + (H z -H x ) 2 ] 

Z£ 2 (2m-\) 2 , „ , o 15-33m-mn-n, „ 6(i-™) 

= — r TTT-intt + C! )-*-— -(n£t + Cl ) 

(m + l) 2 4(9m - 3 + mn + n) 

3c 3 £(l - 2ra) , . , 3 , (m + 1) 2 . „ , r igo-mj 2 i 



m+1 v ; £ 2 (9m-3 + mn + n) 2 

+ 



2C3(m + 1} („ft + cO-t™- 1 ! + c 2 (^ + Cl )"* • (37) 



£(9m - 3 + mn + n) 

V = ABC = (n£t + ci)» , (38) 

2ct 2 _2(2m-l) 2 15-33m-mra-n H 1 _^i +2 



3ff 2 (m + 1) 2 6l 2 (9m-3 + mn + n) 



(n^ + ci) "("+!)" 



2c 3 (l-2m), n-s,-, 2(m + l) 2 . r i 2 (i-m) ^ 



We observe that all the parameters diverge at t = —ci/n£. Therefore, the model has a singularity at 
t = —Ci/n£, which can be shifted to t = by choosing c\ = 0. This singularity is of Point Type as all the 
scale factors vanish at t = —c\/n£. The cosmological evolution of Bianchi-II space-time is expansionary since 
all the scale factors monotonically increase with time (see, Fig.l). So, the Universe starts expanding with a 
big bang singularity in the derived model. The parameters p, p, p p and A start off with extremely large values. 
In particular, the large values of p p and A in the beginning suggest that strings dominate the early Universe. 
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Figure 2: Particle energy density (p p ) and string tension density (|A|) versus time (t). 
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Figure 3: Plot of j& versus time (t). 
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Figure 4: Plot of displacement vector (f3) versus time (t). 
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Figure 5: Plot of density parameter (Q) versus time (t). 



For sufficiently large times, p p and A become negligible. Therefore, the strings disappear from the Universe for 
larger times. 

From Fig. 1, it is observed that during cosmic expansion, C(t) > B(t) > A(t). Thus the derived model can 
be utilised to describe the anisotroy of early Universe. In particular for m = 0.5, at late time, the directional 
scale factors vary as 

A(t) « (nlt + ci)n, B(t) « (n£t + ci)», C(t) « (n£t + ci)» . 

Therefore, isotropy is achieved in the derived model for m = 0.5. 
The value of DP (q) is found to be 

q = -l + ^ = n-l, (40) 

which is a constant. A positive sign of q, i.e., n > 1 corresponds to the standard decelerating model whereas 
the negative sign of q, i.e., < n < 1 indicates acceleration. The expansion of the Universe at a constant rate 
corresponds to n = 1, i.e., q = 0. Also, recent observations of SN la [43]- [47] reveal that the present Universe is 
accelerating and value of DP lies somewhere in the range — 1 < q < 0. It follows that in the derived model, one 
can choose the values of DP consistent with the observations. For n < 1, the model is accelerating whereas for 
n > 1 it goes to decelerating phase. In what follows, we compare the two modes of evolution through graphical 
analysis of various parameters. 

Fig. 2 depicts the variation of p p and |A| versus time in the decelerating and accelerating modes of Universe. 
It is observed that in the early Universe p p > |A|, and for sufficiently large time,p p and |A| tend to zero. Therefore, 
the strings disappear from the Universe at late time (i. e. present epoch). According to Ref. |17j . since there is 
no direct evidence of strings in the present-day Universe, we are in general, interested in constructing models of 
a Universe that evolves purely from the era dominated by either geometric string or massive strings and ends 
up in a particle dominated era with or without remnants of strings. Therefore, the above model describes the 
evolution of the Universe consistent with the present-day observations, has been graphed in Fig. 3. It is 
evident that is decreasing function of time and pjr > 1 in the early Universe. This shows that massive strings 
dominate the early Universe. But the string tension density drops to zero for sufficiently large times. That is 
why the strings are not observed in current astronomical observations. 

From equation (I28[) . it is clear that displacement vector (/?) is decreasing function of time and approaches 
to small positive value. This behaviour of j3 is clearly shown in Fig. 4, in accelerating and decelerating phase 
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of the Universe. From equation it is observed that for n < 3 and m = 0.5, the density parameter (f2) 

approaches to 1 at late time. Thus the derived model predicts a flat Universe at the present epoch. Fig. 5 
depicts the variation of density parameter versus time during the evolution of Universe. 



3.2 Lyra's Cosmology of String fluid with Exponential-law 

Solving the equations (fl7|) . (fT2|) and (Tl9|) . we obtain the metric functions as 



= C4C^ +1 exp 



3m£ c 
-t + 



6(m-l) 

2 



m+1 18^ 2 (3m - l)(m - 1) 



(m + 1) 2 e ££i21_i) t _ _C3_ e _ 3ft 



3£ci 



-4 L 2 



exp 



C(i) = 



Equation (|T6|) leads to 



Integrating equation (f35|) . we obtain 



where fco is an integrating constant. 



3(1 -m)£ 



6(m-l) 

2 



(m + l) 2 ee(.m-i) t C3 



m+1 18£ 2 (3m- l)(m- 1) 

/3 



e "+ 1 



-3£t 



3^3 







-31 



/3 = fco e 



-3£t 



(41) 
(42) 
(43) 

(44) 
(45) 



The expressions for the isotropic pressure, the proper energy density, the string tension and the particle 
density for the derived model are obtained as 



9r(m 2 - m + 1) c 2 m+1 (to + 1) m m - 1 )_ t 2c 3 c 2 m+1 (to + 1 



(to + 1) 2 4(3to-1) 



e "+ 1 



3^(3to - 1) 



-e "+ 1 



12(m-l) 

3c 3 (l - 2m) e _3fl _ c 2 m+1 ( m + l) 2 ^ i2<jm-i) t _ c| R _ 6tt _ 3 fc2 ^ 



4{m + l) 9£ 2 (3to-1) 2 



'-2 



0' 



(46) 



9 = 



6(m-l) 3(m-3) 

9£ 2 to(2 - to) ( c 2 m+1 (5 - 11m) w(m-i) f 2c3C 2 m+1 (m + 1) 3i(m-3) f 



(m + 1) 2 4(3m-l) 



3f (3m - 1) 



12(m-l) 



e + 4 fc e 



3c 3 (l - 2m) 3£t c 2 m+1 (m + l) 2 p 1 4 p -m , 3 z„2, .-»/ 

c|(m+l) 9€ 2 (3m-l) 2 

_ 9^ 2 (2m - 1) 

A — 

m+1 



(47) 
(48) 



9v 



9r(l + m-3m 2 ) t c 2 m+1 (13 — 25m) 2c 3 c 2 m+1 (m + 1) 3*(m-3) t 



+ 



(m + l) 2 4(3m - 1) 

12(m-l) 

3c 3 (l-2m) 3£t c 2 m+1 (m+ l) 2 i2^(^-d t c 



-e ™+! 



3£(3m - 1) 



cf(m+l) 9^ 2 (3m-l) 2 
The critical energy density (p c ) and density parameter (fi) are given by 



e "-" + -A: n e 
Co 4 



2 -3tt 



(49) 



(50) 
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3m(2 — m) c 2 m+1 (5 — 11m) 6*(m-i) t 2c 3 c 2 m+1 (m + 1) 3<(m-3) f 
(m+1) 2 + 12^ 2 (3m-l) 6 ™ +1 9^ 2 (3m- 1) 6 

12(m-l) 

3c 3 (l -2m) 3^ _ c 2 m+1 (m + 1) 2 i2^-D t _ _c|_ Mt J_ , 2 3 ft rM v 

+ 3£ 2 c3(m + l) 27^ 4 (3m-l) 2 3^ 4£ 2 1 j 

The average Hubble's parameter, expansion scalar, shear scalar spatial volume (V), and anisotropy parameter 
(A) of the model are, respectively given by 

H = £. (52) 
6 = U, (53) 

6(m-l) 3(m-3) 

2 3f 2 (2m-l) 2 c 2 m+1 (5 -11m) 6<(m-i) t 2c 3 c 2 m+1 (m + 1) 3i t m-3) t 

(m + 1) 2 4(3m- 1) 6 + 3^(3m - 1) 6 

12(m-l) 

3c 3 (l-2m) 3ft c 2 "' +1 (m + 1) 2 i2<(m-i) t c§ Mf 
c|(m + l) + 9^ 2 (3m-l) 2 cf ' 1 ' 

V = 4e m , (55) 



T 2(2m-l) 2 c, m+1 (5 -11m) 4c 3 c, m+1 (m + 1) gg(m-3) , 

^4 = — 1 £ £ i i e m + 1 1 -I i L e m + 1 1 

(m + 1) 2 6£ 2 (3m-l) 9F(3m-l) 

12(m-l) 

2c 3 (l - 2m) ^_ 3 ft 2c 2 m+1 (m + l) 2 ig^ij t 2c| 6ft 



£ 2 c^(m + l) 27^ 4 (3m-l) 2 ii 2 
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The value of DP (q) is found to be 

q = -l. (57) 

Recent observations of SN la [43] — [47] suggest that the Universe is accelerating in its present state of 
evolution. It is believed that the way Universe is accelerating presently; it will expand at the fastest possible 
rate in future and forever. For n = 0, we get q = — 1 ; incidentally this value of DP leads to dH/dt = 0, which 
implies the greatest value of Hubble's parameter and the fastest rate of expansion of the Universe. Therefore, 
the derived model can be utilized to describe the dynamics of the late time evolution of the actual Universe. 
From equation (|5"Tj) . it is clear that the density parameter (f2) approaches to 1 at late time. Thus the derived 
model predicts a flat Universe. 



4 Concluding Remarks 

In this paper, we have studied Bianchi-II string cosmological models in normal gauge for Lyra's manifold with 
constant deceleration parameter, considering two cases, (3.1) and (3.2). for n^O and n = respectively. It is 
observed that in both cases, the mean anisotropy parameter (A) drops to zero for m = 0.5 for sufficiently large 
times. Thus for m = 0.5, isotropy is achieved in derived models at late time. The main features of the work are 
as follows: 



• The models are based on exact solutions of modified Einstein's field equations in normal gauge for Lyra's 
manifold for the anisotropic Bianchi-II space-time filled with massive strings. 

• The singular model (n 7^ 0) seems to describe the dynamics of Universe from big bang to the present epoch 
while the non-singular model (n = 0) seems reasonable to project dynamics of future Universe. 

• The Universe acquires flatness at late time (see Fig. 5) which is consistent with the prediction of current 
observations. 

• The displacement vector (/3) is a decreasing function of time and it approaches to small positive value at 
late time. This matches with the nature of cosmological constant (A). 
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• The string tension density vanishes for sufficiently large times (see Fig. 2). Thus the strings disappear 
from Universe at late time. That is why, the strings are not observable in the present Universe. 

• The age of Universe in singular model is given by 

which differ from present estimation i. e. T = Hq w lAGyr. But if we take, n = 1 and c\ = 0, the 
model is in good agreement with present age of Universe. 
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